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Recent theoretical works have demonstrated the realization of fractional quantum anomalous Hall 

states (also called fractional Chern insulators) in topological flat band lattice models without an 

. . . external magnetic field. Such newly proposed lattice systems with topological flat bands play a 

CO , vital role to obtain a large class of fractional topological phases. Here we report the exact numerical 

studies of edge excitations for such systems in a disk geometry loaded with hard-core bosons, which 

^3 ' will serve as a more viable experimental probe for such new fractional topological phases. We find 

Cn \ convincing numerical evidence of a series of edge excitations characterized by the chiral Luttinger 

$_i( ■ liquid theory for the bosonic fractional Chern insulators in both the honeycomb disk Haldane model 

Oh' and the kagome-lattice disk model. We further verify these current-carrying chiral edge states by 

•^r \ inserting a central flux to test their compressibility. 

^ ■ PACS numbers: 73.43.Cd, 05.30.Jp, 71.10.Fd, 37.10.Jk 
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Introduction. — One of the most essential and fasci- 
nating topics in condensed matter physics is to explore 
and classify the various states of matter, among which 
the integer quantum Hall effect (IQHE) [1| and the frac- 
tional quantum Hall effect (FQHE) y have long been 
the major focus. It is well known that, at fractional 
fillings, the FQHE will emerge when interacting parti- 
cles move in Landau levels (LLs) caused by an external 
uniform magnetic field. On the other hand, great in- 
terest has been aroused to realize both the IQHE and 
FQHE in lattice models in the absence of an external 
magnetic field. An initial theoretical attempt was made 
by Haldane 0], who proposed a prototype lattice model 
to achieve the IQHE by introducing two non-trivial topo- 
logical bands with Chern numbers C = ±1 [J]. Haldane's 
model demonstrates that the IQHE can also be attained 
without LLs, i.e. defines the quantum anomalous Hall 
(QAH) states. The lattice version of the FQHE without 
LLs, however, comes much latter because of its intriguing 
strongly correlated nature. 

The recent proposal of topological flat bands (TFBs) 
fulfils the basic ingredients to explore the intriguing frac- 
tionalization phenomenon without LLs. TFBs [5|-|7| be- 
long to a class of band structures with at least one nearly 
flat band with non-zero Chern number. This may be 
viewed as the lattice counterpart of the continuum LLs. 
Several systematic numerical works have been done to 
explore the correlation phenomenon within TFB mod- 
els, and AbeUan ^^ and non-Abclian [uHlll FQHE 
states without LLs have already been well established 
numerically. This intriguing fractionalization effect in 
TFBs without an external magnetic field, defines a new 
class of fractional topological phases, i.e. fractional quan- 
tum anomalous Hall (FQAH) states, now also known 
as fractional Chern insulators (FCIs). Some new ap- 
proaches, e.g. the Wannier-basis model wave functions 



and pseudo-potentials [14|, the projected density opera- 



tor algebra 15[ , the parton wave- function constructions 
jla |. and also the adiabatic continuity paths [l7|, have 
been proposed to further understand these FCI/FQAH 
states in TFBs. There are various other proposals of TFB 
models and material realization schemes |18l434| . Very re- 
cent systematic numerical studies have also found exotic 
FCI/FQAH states in TFBs with higher Chern numbers 
(C > 2) which do not have the direct continuum anal- 
ogy in LLs [35l - [37t . and also the hierarchy FCI/FQAH 
states 1381. 



Although previous studies have firmly established var- 
ious properties of a large class of FCI/FQAH phases, we 
are here concerned with the less studied edge excitations 
since they in principle should open another window to 
reveal the bulk topological order [39|. Edge excitations 
might also provide a more viable experimental probe (40[ 
especially considering about the possible future realiza- 
tions of FCI/FQAH phases in optical lattice systems [30| . 
A recent work has studied edge excitations of the bosonic 
FQHE in LLs generated by an artificial uniform flux in 
optical lattices [41|. In the present work, we investigate 
edge excitations of the bosonic FCI/FQAH phases 19| of 
TFB models in a disk geometry. Through extensive sys- 
tematic numerical exact diagonalization (ED) studies, we 
demonstrate clear evidence of edge excitation spectra of 
bosonic FCI/FQAH phases. These edge excitation spec- 
tra are quite coincident with the chiral Luttinger liquid 
theory [39| . To further check the compressibility of these 
edge states, we insert a central flux into the disk sys- 
tem and indeed verify these current-carrying chiral edge 
states upon tuning the flux strength. 



Model Hamiltonians. — We first look into the Haldane 
model [3[ on a honeycomb-lattice disk, which is loaded 




FIG. 1: (color online), (a) The honeycomb-lattice disk and 
(b) the kagome-lattice disk. Both disk geometries satisfy the 
Ce rotational symmetry. Various disk sizes are indicated by 
the circles and the labelled numbers. 



with interacting hard-core bosons [9[ : 
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where 6j creates a hard-core boson at site r, (...), {{■■■)) 
and (((...))) denote the NN, the NNN and the next- next- 
nearest-neighbor (NNNN) pairs of sites, respectively. We 
adopt the previous parameters [9| to achieve a lowest 
TFB with a flatness ratio of about 50 (the ratio of the 
band gap over bandwidth): t = I, t' = 0.60, t" = -0.58 
and 4> = OAtt. 

We also consider a kagome lattice TFB model [5|, |25[ 
which is also loaded with hard-core bosons, and the 
model Hamiltonian has the form: 
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We choose t he p revious parameters for this kagome- 
lattice model ^: t = 1, t' = -0.19, (j> = 0.227r, which 
leads to a lowest TFB with the flatness ratio of about 20. 

Convincing numerical evidence of the 1/2 bosonic 
FCI/FQAH states at the 1/2 hard-core boson filling of a 
lowest TFB on a torus geometry has been well established 
previously [9[, with the formation of a quasi-degenerate 
ground-state (GS) manifold, characteristic GS momen- 
tums, a robust bulk excitation spectrum gap, as well as 
a fractional quantized Chern number for each GS. The 
disk geometries for both lattice models are illustrated in 
Fig. [U which shows the Cg rotational symmetry. An 
additional trap potential is required on the finite-size 
disk systems to confine the FCI/FQAH droplet, outside 
which the edge modes are able to propagate around the 
disk. Here we choose the conventional harmonic trap, of 
the form V = Vtrap J2r kP'^r with Vtrap as the potential 
strength, |r| as the radius from the disk center, and n^ 
as the boson number operator. 
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FIG. 2: (color online). Edge excitations on a honeycomb disk 
with Na = 96 sites and trap potential Vtrap = 0.005. Numbers 
labelled upon low energy levels in each sector show the quasi- 
degeneracy of low edge excitations, which is also the sequence 
derived from the chiral Luttinger liquid theory. 



Edge excitations. — We have first considered the hon- 
eycomb disk model, and have observed clear edge excita- 
tion spectra for a relative broad range of the trap poten- 
tial strength T4rap and various honeycomb disk sizes with 
Ng = 24, 54, 96 sites. Some representative resuhs are 
shown in Fig. [5] As the system owns the Cg rotational 
symmetry, each energy states can be classified with a 
quantum number of angular momentum L = 0, 1, 2, 3, ..., 
mod 6. These edge-excitation quasi-degeneracies of a fi- 
nite system with various boson numbers turn out to be a 
good match with that predicted by the chiral Luttinger 
liquid theory at least for the lowest six sectors. 
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TABLE I: Occupation of chiral edge bosons in fc = 
1,2,3,4,5,6, ... orbitals: {ni, n2,n3, 724,715, ng, ...}. d is the 
degeneracy of the edge excitations in a sector. AL — L — 
Las = Ylk ^'^'^ ' which is the shifted total angular momentum 
relative to the ground state angular momentum. The right 
column shows the degeneracy in systems with finite boson 
numbers, e.g., for the "4b" column {Nb — 4), the predicted 
degeneracy is 1, 1, 2, 3, 5, 6, 9, ... in AL — 0, 1, 2, 3, 4, 5, 6, ... an- 
gular momentum sectors. The partial degeneracy sequences 
are mostly observed (except the underlined 9, 10 and ]_1) in 
our ED results. 



According to the hydrodynamical approach [39[, low 
energy edge excitations of a Laughlin-state FQHE 
droplet at filling v = l/m are generated by the Kac- 
Moody algebra, and form a chiral Luttinger liquid with 
the effective Hamiltonian: 
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with the angular momentum k along the edge and 
the Fourier-transformed one-dimensional density opera- 



son/phonon operators. Based on this theory, chiral edge 
bosons occupy in A: = 1,2,3,4,5,6,... (unit 2tt/M and 
edge length M) angular momentum orbitals with the oc- 
cupation numbers (711,712,713,714,715,716, ...}. We denote 
AL = L — Lgs a-s the shifted total angular momentum 
relative to the GS angular momentum. For a given AL = 
X^fc^JT-fc, it is easy to count the edge-state degeneracies 
by numerating the allowed occupation configurations. In 
the thermodynamics limit with infinite boson numbers, 
the degeneracy sequence should be 1,1,2,3,5,7,11,15,22,... 
For our systems with only small boson number iV^'s, edge 
excitations made up of ^j, Uk < N}, single modes are 
expected. For three bosons (iVb = 3), after discard- 
ing some configurations in the middle column of Table 
Uwith ^^, Uk > Nb = 3, the predicted degeneracy should 
be 1,1,2,3,4,5,7,..., which is in good accordance with the 
observed quasi-degeneracy in our ED results [Fig. [2jb)]. 
For four bosons {Nb = 4) , after discarding some configu- 
rations with J2k nk > Nb = 4, the predicted degeneracy 
should be 1,1,2,3,5,6,..., which is also in good accordance 
with the observed quasi-degeneracy in our ED results 
[Fig. [UJc)] . The right column in Table U list the par- 
tial degeneracy sequences which are mostly observed in 
our ED results. 

Now we turn to investigate the kagome disk model. We 
studied the finite system also in several different sizes, 
with Ns = 12,30,42,72 sites. Edge spectra from a fi- 
nite kagome disk with Ng = 72 sites are shown in Fig. 
[31 which shows even clearer edge excitation spectra than 
that obtained from Ns = 96 honeycomb disk. It is quite 
exciting to observe that these edge excitations are inde- 
pendent of lattice geometry even in such small systems. 
For an example with six bosons (Nb = 6), after discard- 
ing some states with J^k ^k > Nb — 6, the predicted 
degeneracy should be 1,1,2,3,5,7,... which is also in ex- 
act accordance with the observed quasi-degeneracy in our 
ED results [Fig. [Hd)]. It is also observed that with the 
larger disk sizes and smaller boson numbers, the num- 
ber of matched momentum sectors increases, and thus 
displays less severe finite size effects. 

Central flux. — As a next step, we insert a flux into 
the center of disk geometry and tune the flux strength 
to test the compressibility of edge excitations. Consider 
the kagome disk with Ng = 42 sites for an example; the 
trap potential will always be fixed as Vtrap = 0.005. For 
systems with Nb = S bosons, low energy levels En in 
k = 0,3 sectors (or A: = 1, 4 and fc = 2, 5 sectors) evolve 
into each other, as is shown in Fig. IDJa). These energy 
spectra £"„ finally return to themselves after two peri- 
ods of the central flux, i.e. En{^ + 47r) = En{^). For 
systems with Nb = 4 bosons, low energy levels En in 
fc = 0,2,4 sectors (or fc = 1,3,5 sectors) evolve into 
each other, as is shown in Fig. Wlb)- These states 
finally return to themselves after three periods of the 
central flux, i.e. i?„(<f> -I- 67r) = _£„(<!>). For our disk 
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FIG. 3: (color online). Edge excitations on a kagome disk 
with A'^s ~ 72 sites and trap potential Vtrap = 0.005. Numbers 
labelled upon low energy levels in each sector show the quasi- 
degeneracy of low edge excitations, which is also the sequence 
derived from the chiral Luttinger liquid theory. 



systems with the Cq rotational symmetry, vifhen boson 
numbers satisfy -/V;,/6 = p/q with p and q as coprime 
integers, there is a generic periodicity of energy spec- 
tra: £'„($ + 2g7r) — £'„($). It is obvious to observe 
the current-carrying nature, i.e. the compressibility of 
edge states since all of them evolve into the higher en- 
ergy spectra without any gap. These low edge excita- 
tions are in principal all degenerate in the thermody- 
namic limit, as is explained in the microscopic theory of 
FQHE droplet. When the boson numbers increases in 
the disk system, mixture of bulk topological insulating 
states become apparent. For the system with Ng = 42 
sites and Nf, = 6 bosons, we thus get an incompressible 
ground state, which is always gapped from other excited 
states as shown in Fig. HJ^c). 

Summary and discussion. — We considered two repre- 
sentative TFB lattice models in disk geometry to explore 
edge excitations of FCI/FQAH states. With a confin- 
ing harmonic trap potential, very clear edge excitation 
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FIG. 4: (color online). Evolution of low energy states with 
the varying central flux of the kagome disk with Ns = 42 sites, 
trap potential Vtrap ~ 0.005, and various boson numbers Nt^s. 
(a) For A^i, = 3, low energy states in fc = 0, 3 sectors (red, 
green respectively) evolve into each other, (b) For A'^;, = 4, 
low energy states in fc = 0, 2, 4 sectors (red, green, blue re- 
spectively) evolve into each other, (c) For A^j, = 6, low energy 
states in all six sectors are displayed, and the ground state is 
gapped from other excited states, showing the incompressible 
nature of an insulating state. 



spectra are observed with the quasi-degeneracy count- 
ing rule satisfying the chiral Luttinger liquid theory at 
least for the lowest six sectors. By inserting and tuning 
a central flux to both systems, we are able to examine 
the compressibility of these current-carrying chiral edge 



states. Our work here focuses on the bosonic FCI/FQAH 
states in TFBs with the Chern number C = 1 , while it is 
very interesting to explore the exotic ed ge e xcitations in 
high-Chern-number FCI/FQAH phases 13514371 and also 
those in the hierarchy FCI/FQAH phases [38[ in future 
works. We also expect that edge excitation spectra pro- 
vide another window to reveal the bulk topological order, 
as well as a more viable experimental probe to future 
FCI/FQAH experimental systems. 
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